The existence, nonexistence and uniqueness of 
global positive coexistence of a nonlinear 
elliptic biological interacting model 

Joon Hyuk Kang, Yun Myung Oh 
Department of Mathematics, Andrews University 
Berrien Springs, MI 49104, U.S.A. 
(kang@andrews.edu, ohy@andrews.edu) 

Abstract 

The purpose of this paper is to give a sufficient condition for the 
existence, nonexistence and uniqueness of coexistence of positive so- 
lutions to a rather general type of elliptic competition system of the 
Dirichlet problem on the bounded domain £1 in R n . The techniques 
used in this paper are upper-lower solutions, maximum principles and 
spectrum estimates. The arguments also rely on some detailed prop- 
erties for the solution of logistic equations. This result yields an alge- 
braically computable criterion for the positive coexistence of compet- 
ing species of animals in many biological models. 
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1 Introduction 

The coexistence of steady states of competition interacting models with dif- 
fusion has been an object of intensive study in recent years. See, for example, 
lists of references in [1], [3], [4], [5], [9], [10], [14], [15], [14], [15]. The 
most general type of parabolic competition interacting system is 

u t = Au + ug(u, v), 
v t = Av + vh(u, v), 

where A is the Laplacian and u, v represent the densities of two competing 
species of animals. The terms Am and Av model dispersal by means of simple 
diffusion. We assume here that the C 1 functions g and h are relative growth 
rates satisfying the following so-called growth rate conditions: 

(Gi) g u (u, v) < 0, g v (u, v) < 0, h u (u, v) < 0, h v (u, v) < 0, 

(G2) There exist constants Co > 0, c\ > such that g(u, 0) < for u > Co 

and h(0, v ) < for u > C\. 

The hypothesis [G\) characterizes how the two species u and v interact 
with each other in terms of their relative growth rates. It is well known that 
condition (G2) exhibits the so-called logistic pattern while the constants Co 
and ci are referred to as the carrying capacity. 

The earlier literature on this line focused in the Neumann boundary value 
problem: 

u t = Au + ug(u, v), 

v t = Av + v h(u, v), (1) 

t = = t for (*.*) e M x d ^ 

and in its steady state, the elliptic system 

A-u + ug(u, v) = 0, 
< Av + vh(u, v) = 0, (2) 
duj*l = = Me) for x e dtt, 

v an on 1 

where n denotes the unit out-normal along boundary dfl. The Neumann 
boundary conditions = = are interpreted as an assumption that 
both populations are staying inside, that there is no migratory flux across 
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dfl. The goals of investigations along this line include finding out under what 
conditions on the nonlinearities g and h systems (1) and (2) have positive 
solutions u > 0, v > and the possible uniqueness. Most of the work in this 
case were established by P. DeMottoni and F. Rothe in 1979 [7] and P. Brown 
in 1980 [2]. Their work in a large sense completes the avenue of investigation 
in the study of Neumann boundary value problems. Researchers thus have 
since turned their attention to the biologically and physically more important 
case that is the Dirichlet boundary condition: 

Au + ug(u,v) = 0, 
< Av + vh(u,v) = 0, (3) 
(u,v)\ dn = (0,0). 

Biologically, this setting allows migration of these two populations across the 
boundary but they may not stay on dQ, where, for example dfl is a river. 
It was then found that the features known in the Neumann setting are not 
usually shared by those in the Dirichlet setting. The study in the latter 
setting, especially in the case of steady states like system (3), seems to be 
more difficult. 

The goal of this paper is to answer the following questions about positive 
steady state to (3). 

Problem 1 : What are the sufficient conditions for existence of steady state? 
Problem 2 : Is it possible for either one of the species to be extinct? 
Problem 3 : When is the coexistence state unique? 



2 Preliminaries 

In this section we state some preliminary results which will be useful for our 
later arguments. 

Definition 2.1 (Upper and Lower solutions) 

The vector functions (u 1 , ...,u N ), (u 1 , ■■■,u N ) form an upper/lower solution 
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pair for the system 



Au l + g' l (u\...,u N ) = in Q 
u % = on dVL 



if for i 



N 



1 1 ,u l ,u l+1 



Au % + g l (v}, u 
Au l + g^u 1 , u\ u 



i • • • i 

i+l „,N\ 



i • •• i 



U 



and 



> 

in Q for u? <u 3 < u> , j ^ i, 

u l <u % on Q 

u l < < u l on dVt. 



Lemma 2.1 ( flj) 

If g l in the Definition 2. 1 are in C 1 and the system admits an upper/lower 
solution pair (u 1 , ...,u N ), (u 1 , ...,u N ), then there is a solution of the system 
in 2.1 with u l < u l < u l in Cl. If 

Ae + g\u 1 ,...,u N ) ^ 0, 
Am* +g i (u 1 ,...,u N ) ^0 

in Q for i — 1, then u % < u % < u l in Q. 

Lemma 2.2 (The first eigenvalue) ([6]) 

j —Au + q(x)u = \u in f2, 
| u\ dn = 0, 



(4) 



where q(x) is a smooth function from Q to R and Q is a bounded domain in 
R n . 

(A) The first eigenvalue Xi(q) of (4), denoted by simply Ai when q = 0, is 
simple with a positive eigenf unction. 

(B) If q±(x) < q2(x) for all x G Vt, then Ai(gi) < Ai(g 2 )- 

(C) (Variational Characterization of the first eigenvalue) 



\\{q) = min 



Jg(|V0| 2 + # 2 )^ 
In <P 2 dx 
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We also need some information on the solutions of the following logistic 
equations. 

Lemma 2.3 ([14]) 

{ Au + uf(u) = in Q, 
\ u\ dn = 0,u > 0, 

where f is a decreasing C 1 function such that there exists Co > such that 
f(u)<0 for u > c and Q is a bounded domain in R n . 

(1) If f(0) > Ai, then the above equation has a unique positive solution, where 
Ai is the first eigenvalue of —A with homogeneous boundary condition. We 
denote this unique positive solution as Of. 

(2) If f(0) < X\, then the above equation does not have any positive solution. 

3 Existence, Nonexistence and Uniqueness 

We consider the system (3) with conditions (Gi) and (G 2 ). 

Theorem 3.1 (A) If g(0, c() > Ai and h(c ,0) > X±, then (3) has a solution 
(u, v) with 

ff (.,ci) <u< e g{ . fi) 

Qh(co,-) < V < 9h(0,-)- 

Conversely, any solution (u,v) of (3) with u > 0,v > in Q must satisfy 
these inequalities. 

(B) If 5(0,0) < Ai or h(0,0) < X\, then (3) does not have any positive 
solution. 

Proof. 

(A) Let u = ^ s (.,o)) v = 0h(o,-)- Then by the monotonicity of g, 

Au + ug(u, v) 
= Au + u(g(u,0) - g(u,0) +g(u,v)) 
= u(g(u,v) -g(u,0)) < 0. 

Similarly 

Av + vh(u, v) < 0. 
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So, (u,v) is an upper solution to (3). 

Let u = 9 g (., ci ) and v = 0h(c o ,-)- Then by the Maximum Principles, we obtain 



u < Og(.fl) < c , 
v < 0ft(o,-) < ci. 



By the monotonicity of <7, 



Am + ug(u, v) 
= Au + u(g(u,ci} - g(u,ci} + g(u,v)) 
= u(g(u,v) -g(u,d)) > 0. 

Similarly, 

+ vh(u, v) > 0. 

Therefore, (u,v) is a lower solution to (3). Furthermore, u < u,v < v in f2 
and M = -u = v = v = 0on cftl 
So, (3) has a solution (w, v) with 

g( . ;Cl) < u < fl( . jO ), 
^(c ,-) < v < ^(o,-)- 

Suppose (m, is a coexistence state for (3). Then since 

Au + ug(u, 0) 
> Au + ug(u,v) = 0, 



u is a lower solution of 



AZ + Zg(Z,0) = in Q, 
Z = on 



(5) 



But, since any constant larger than c is an upper solution of (5), we have 

u < 9 g{ . :0) . (6) 

Similarly, we have 

v < O h{0 .). (7) 
Since v < 9h(o,-) < c i, by the monotonicity of g 

Au + ug(u, ci) 
< Am + «</(«, ?;) = 0. 
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Therefore, u is an upper solution of 



AZ + Zg(Z, Cl ) = in Q, 

Z = on da { ' 

If e > is so small that g(e(j)i, c±) > Ai on Q, where 0i is the first eigenvector 
of —A with homogeneous boundary condition, then since 

Ae0i + e0i#(e0i,ci) 
= e(A0i + 0i5(e0i,ci)) 
> e(A0! + Ai0i) = 0, 

e0i is a lower solution of (8). So, we have 

^(•.d) < u. (9) 



Similarly, we have 

By (6), (7), (9) and (10), 



0fc(co,O < ( 10 ) 



ff (.,ci) < M < ^(-,0), 

^(co,0 < v < 0/i(o,o- 



(S) Assume #(0,0) < A x . The other cases are proved similarly. Suppose 
(u,v) is a positive solution to (3). Then since 

Au + ug(u, 0) 
= Au + u(g(u,v) -g(u,v) + g(u,0)) 
= u(g(u,0) -g(u,v)) > 0, 

w is a lower solution to 

AM + «g(«,0) = in f2, , s 

u = o on aa ^ j 

Any constant larger than Co is an upper solution to (11). Hence, (11) has a 
positive solution -u with u < u . This contradicts to the Lemma 2.3 which 
says there is no positive solution of (11) if g(0,0) < X±. 



7 



Theorem 3.2 lfg(0,a) > Ai,/i(c ,0) > Ai and 

4mf(-^)inf(-^ i ) > ^(sup^) 2 + ^(sup^ 1 ) 2 

+2(sup^)(sup^), 

then (3) has a unique positive solution. 



Proof. Suppose (ui,vi) and (u 2 ,v 2 ) are positive solutions to (3). Let 
p = u\ — U2 and q — Vi — v 2 . Then 

Ap + p</(ui,ui) 
= Atti - Au 2 + (tii - u 2 )g(u 1 , uj) 
= -Am 2 - u 2 g(u!, v x ) 
= -Au 2 - u 2 (g(u 2 , v 2 ) - g(u 2 , v 2 ) + 

= -U2(p(ui,Ul) - g(u 2 ,V 2 )) 

= -u 2 (g(u!, ui) - c/(w 2 , ui) + 0(ii2, «i) - y(«2, 

But, by the Mean Value Theorem, there is z depending on iti, w 2 such that 

/ \ / \ dg(x,V!) 
g{u u Vi) - g{u 2 , vx) = — — — p. 



Hence, 



UQ\X V\j 

Ap + p</(ui,ui) = -u 2 [ — — p + g(u 2 ,vi) -g(u 2 ,v 2 )\. 



i.e., 



Ap + c/(^i, vi)p + u 2 p^§^- - u 2 (g(u 2 , v 2 ) - g(u 2 , v x )) = 0. (12) 



du 

The same argument shows that 

Aq + h(u 2 , v 2 )q + V!q dh{ ^' x) - Ui(fr(u 2) «i) - /i(ui, v i)) = °> ( 13 ) 

where x depends on i>i , u 2 by the Mean Value Theorem. 

Since \i(—g(ui,Vi)) = 0, by the Variational Characterization of the first 

eigenvalue, 

f sl Z(-AZ-g(u 1 ,v 1 )Z)dx>0 (14) 
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for any Z G C 2 (Q) and Z\qq = 0. The same argument shows that 

J Q W(-AW -h(u 2 ,v 2 )W)dx > (15) 

for any W G C 2 (fi) and W\ an = 0. 
From (12) and (13), we get 

-pAp - g(m, vi)p 2 - d9i y i] u 2 p 2 + u 2 p(g(u 2 , v 2 ) - g(u 2 , v ± )) = 0, 
-qAq - h(u 2 , v 2 )q 2 - &h ^l' x) Viq 2 + V!q(h(u 2 , v ± ) - h(u 1} vj) = 0. 

Hence from (14) and (15), 

In(- d9i du l) u ^P 2 + u M9(u 2 , v 2 ) - g(u 2 , v^) + viq(h(u 2 , - h(u u v^) 
-^ff^ Vl q 2 )dx < 0. 

By the Mean Value Theorem, for each x G Q, there exist y, y such that 

g(u 2 ,v 2 )-g(u 2 ,v 1 ) = d -^(-q), 
h(u 2 ,v 1 )-h(u 1 ,v 1 ) = 9 -^pl(-p), 

which implies that 



/ 

Jn 



dg(x,vi) 2 dg{u 2 ,y) dh{y,vi) dh(u 2 ,x) 2 

w 2 p - («2 — « h *>i — 5 )p<? 5 wi<? fife < 0. 



In du "' dv du ' dv 

Therefore, we find 

_ _ n - f dg(x,V!) 2 . dg(u 2 ,y) dh(y,Vi). 

p = q = if u 2 ( - (u 2 + vi )(r] 

du dv du 

V\rj 2 is positive definite 



dv 

for each x G Vt. 



This is the case if 



u^M^Mf + ^ ( ^i) )2 + 2u 2 v l d ^ d f )dh f: i) 
-4 d9 ^ Vl)dh{ ; 2 / ) u 2 v 1 < for each xeQ. 

j e /^ 9g{x,v{) dh(u 2 ,x) y u 2 / dg(u 2 ,y) \ 2 , t>i / 1) \2 

+2 a 9 (u 2 ,y) dhfs.v!) f oreachxG a 
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But, from the inequality in (A) and the hypothesis in the theorem, 

U2 ( dg(v,2,y) \2 , 2 9 9(^2,y) dh(y,vi) , Vl / dh(y,vi) \2 

< ^(sup^) 2 + ^(sup^ 
+2sup(^)sup ' 



(16) 



ov ' L v au ' 

< 4in f(_M^)) in f(_^)) 

dg(u 2 ,x) 
— 9m dv ' 

We can also extend the results to the case when there are multiple species 
competing in the same environment. 

Consider the interacting model 

Aui + u i g i (u i ,u 2 , ...,Ui,u i+1 , ...,u N ) = in Q 

Ui = on <9f2 

for i = 1,...,JV. 

Again, we assume here that the C 1 functions ^ for i — 1, TV are relative 
growth rates satisfying the following growth rate conditions: 

(Ml) g£<0fori,j = l,2,..,iV, 

(M2) There exist constants Ci > 0, c 2 > 0, cjy > such that 
^(0, 0, Ui, 0, 0) < for Ui > d. 

Again, (Ml) characterizes how the N species ui, u 2 , ...,un interact with each 
other in terms of their relative growth rates and (M2) is the logistic pattern 
with carrying capacity constants c±, c 2 , c^. 

The followings are the main results. The proofs are similar to those with 2 
competing species, and so we just sketch it without the details. 

Theorem 3.3 (A) If g^d, c 2 , 0, c i+1 , c N ) > Ai fori = 1,...,N, 
then (16) has a solution (ui, u N ) with 

^p i (ci,...,c i _i,-,c i+ i,...,cjv) < Ui < # 9i (o,...,0,-,0,...,0) 

fori = l,..,N. 

Conversely, any solution (ui,...,un) of (16) with U; L > in Q must satisfy 
these inequalities. 

(B) If <7j(0, 0) < Ai for some i = 1,...,N, then (16) does not have any 
positive solution. 
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Proof. 

{A) Let Ui = gi( o,...,(VA...,o) and Ui = ^(ci,.., Ci _ 1 ,., Ci+1) ..,c JV ) for z = l,...,iV. 
Then by the Maximum Principles and the monotonicity of (u±, u^, u^) 
and (ui, ...,Ui, un) are upper and lower solutions to (16), respectively. 
Furthermore, for i — 1, N, < u~i in Q, and ui = Ui = on dQ. 
So, (16) has a solution (u±, ...,un) with the desired inequalities 

Qgi{ci,-,Ci-i,-,c i+1 ,...,c N ) < Ui < # Si (0,...,0,-,0,...,0) 

for z = 1,..,7V. 

Suppose (ui, ...,un) is a coexistence state for (16). Then by the direct com- 
putation using the monotonicity of g iy we know that Ui is a lower solution 
of 

AZ + Z 9i (0,...,0,Z,0,...,0) = in Q, 

z = o on on { > 

for % = 1, N 

But, since any constant larger than q is an upper solution of (17), we have 

Ui < 9i (o,... ) o,-,o,...,o) (18) 

for i = 1, N. 

Since itj < 9i (o,...,o,.,o,...,o) — c «> by t ne monotonicity of g iy we can derive that 
Ui is an upper solution of 

AZ + Zg i (c 1 ,...,c i - 1 ,Z,c i+1 ,...,c N ) = in Q, , . 

Z = on 50 1 J 

for z = 1, AT. 

If e > is so small that #j(ci, q_i, e0i, c i+ i, c^) > Ai on Cl, where 4>i is 
the first eigenvector of —A with homogeneous boundary condition, then by 
the dirct computation again, we know that e</>i is a lower solution of (19). 
So, we have 

Qg i (ci,...,C i - 1 ,;Ci + i,...,CN) < U i (20) 

for i = 1, N. 
By (18) and (20), 

Ogi(ci,...,Ci-i,;Ci + i,...,c N ) < u i < ^ 9i (0,...,0,-,0,...,0) 

for % = 1,...,N. 
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(B) Without loss of generality, assume g±(0, ...,0) < Ai. 

Suppose (mi, ...,un) is a positive solution to (16). Then by the monotonicity 

of <7j, Ui is a lower solution to 

AZ + Z 9l (Z,0,-,0) = in n, 

z = o on on. { > 

Hence, by the fact that any constant larger than c\ is an upper solution to 
(21), (21) has a positive solution u\ with u x < u\ that contradicts to the 
Lemma 2.3. 

Theorem 3.4 If gi(ci, 0, Cj+i, cn) > Xi and 



W 9l -(0,...,0,-,0,...,0) 

unique coexistence state. 



for i = 1,...,N, where K = sup, 3 9j(0 '"'' ''' 0, '"'" ) , then (16) has a 

,Ji ~ a 9i(c 1 ,...,c i _ 1 ,-,c i+1 ,...,c N ) 



Proof. Suppose (ui,...,un) and (vi,...,vn) are coexistence states of (16) 
and let Wi = Ui — Vi for i = 1, N. Then by the direct computation and the 
Variational Characterization of the first eigenvalu, we obtain 



/ 



N 

J2i v i w i(9i(vi, ...,Vi, v N ) - gi(u u Ui, u N ))]dx < 0. 



n i=i 



By the Mean Value Theorem, there exist f and z' lj such that 

In EilJ ^^ 1 '-'^^ 1 '-^^ -^)^ 



(22) 



If the integrand in the left side of (22) is positive definite, then (22) im- 
plies that Wi = in Q for i = 1, ...,N, which means the uniqueness of the 
coexistence state for (16). But for any e > 0, 

%(M1 " J -;f' , ' i+1 --' %) (-t)> i («; J ) 

< dgiju!,...,^-!^ ,v j+1 ,...,v N ) / ^.U^j 1 
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So, we can see that the integrand is positive definite if for i — 1, ...,7V and 
x G Q, 



dgi(vi,...,Vi-i,t\v i+1 ,...,v N ) / 

dxi ' ' 



N ■£-g i {ui,...,u j - 1 ,z^,v j+1 ,...,v N )(-v i ) e-£rg j {ui,...,Ui- 1 ,zi i ,v i+1 ,...,v N )(-v j ) 

> l^j=l,j^i\ 2e ' 2 / 

or equivalently, 



9gi(^l,...,^-l,f',-»i + l,...,-»jy) 

jV -^f-gi(m,...,^-i,^,^ + i,...^ jv ) 

fslTSj'^li-jWi-lj^Nfi+l,-,^)-^ s 

1 2 ~r 

Since fli ( cl) ... jCi _ 1) . )Ci+lv .. jC j V ) < ^ < ^( ,...,o,-,o,...,o) in ft for i = 1,...,7V, (23) 
will hold if for i = 1,...,7V, 



dffi (i;i,...,^i_l,t%i;i + i,...,Ujv) 

dXi 

v Ar ( sup{ ~9^] , £SU P(-aS7) %(o,..,(v,o, ..,o) 



as, ^ 

g i (c 1 ,...,c i _ 1 ,-,c i+1 ,...,c Ar ) ' 



Let K = sup, , . 93( ° ,,,,,0, ' ,0,,,,,0) . Then (23) holds if 

da, SL su P(-|f Xesup(-|^) 

Choosing e = 1, we have 
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